Interband physics in an ultra-cold Fermi gas in an optical lattice 
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We study a gas of strongly polarized cold fermions in an optical lattice when the excited p-bands 
are populated. We derive the relevant Hamiltonian and discuss the expected phase diagram for 
both repulsive and attractive interactions. In the parameter regime covered here, checkerboard 
anti-ferromagnetic ordering is found to be possible for repulsive interactions while for attractive 
interactions, transitions between different types of paired phases are predicted. 



I. INTRODUCTION 

Experiments with ultra-cold Fermi gases in optical lat- 
tices have opened a way to study experimentally inter- 
acting fermionic systems in a highly tunable environ- 
ment [l|, |3]. Among other things, in such a system one 
can expect a multitude of correlated fermionic phases as 
well as superfluidity. In a two-component fermionic gas 
the atom numbers of different components can be inde- 
pendently controlled and such strongly interacting polar- 
ized fermion gases have been recently studied experimen- 
tally [3, 0,01 • Such studies have revealed, for example, 
intriguing phase separation properties and, depending on 
parameters, appearance of trap physics beyond the local- 
density approximation 6]. 

In this article we address the issue of a strongly polar- 
ized two-component Fermi gas in an optical lattice. Po- 
larized Fermi gases on the lowest band have been studied 
previously, see for example Ref. Q, but here we wish 
to investigate the relevant theory as well as phases as- 
sociated with it, when the majority component fills the 
lowest band and also populates excited bands. To this 
end, we derive the Hamiltonian for this system and apply 
it to discuss the possible phases both for repulsive as well 
as attractive interactions. 

Multiband physics has a long history in the condensed 
matter theory and superconductors with several bands 
have been addressed within the BCS as well as Ginzburg- 
Landau formalism [aQ- However, a system of polarized 
ultracold atoms in optical lattices is different in many im- 
portant respects. First, in the theory of transition met- 
als, for example, the multiband description is motivated 
by the fact that the Fermi surface might pass through 
two (or more) bands. Under such circumstances there 
are processes which transfer electrons between bands as 
well as interactions inside the bands. In the system that 
we consider the bands are well separated and there are 
no processes that transfer atoms between bands. Here 
the interaction between the bands is proportional to the 
product of the densities in different bands and does not 
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involve transfer of atoms between bands 'S|. Second, in 
the present system only interband interactions are rel- 
evant. This makes it much easier to study interband 
effects without being disturbed by large inband effects. 
Third, the anisotropic nature of the p-band tunneling as 
well as the multi-flavor character of the p-band fermions 
gives rise to qualitatively novel possibilities. Fourth, the 
Hamiltonian for the system we discuss can be derived 
from the microscopic theory in a controlled way. 

Higher band physics with Bose-Einstein condensates 
was experimentally studied by Miiller et al. [lO| and inter- 
action induced transitions to higher bands were observed 
by Kohl et al. [1] while Diener and T.-L. Ho [ll| tackled 
this problem theoretically. In this paper, however, the in- 
teraction strengths are much smaller than the bandgaps 
and such interaction induced higher band effects do not 
play an important role. Other interesting possibilities 
also exist. For example, the higher band atoms might 
have stronger nearest neighbor interactions which might 
make supei^ solid phases appear in some parts of the phase 
diagram \1^ . On the other hand when the lattice geom- 
etry is varied, unconventionally ordered quantum phases 
are predicted to be possible [ij, [ij. Some aspects of 
the higher band physics in one dimension were also dis- 
cussed by Kantian et al. [l^ in the context of lattice 
excitons, by Karkkainen et al. 16] for repulsive interac- 
tions with high filling fractions, and by A. F. Ho ,17j], who 
focused on the equal mixture Fermi gas with strong inter- 
actions and with two atoms per site. More recently anti- 
ferromagnetic properties of p-band fermions at half filling 
were studied by Wu and Zhai la\ and p-band fermions in 
a two-dimensional lattice with different lattice geometries 
by Zhao and Liu |19| . 



II. HAMILTONIAN 

Atoms in a cubic optical lattice experience a potential 
K(r) = J2a^a,<Tsin^'^ra/d, where a = {x,y,z}, Va,cr 
are the lattice depths for atoms of type a, and d is the 
lattice constant. Our interest is in a strongly polarized 
two-component gas where the majority component fills 
the lowest band and occupies also part of the first ex- 
cited bands. We label the components by f and |. A 



Hubbard-model Hamiltonian for the fermions is arrived 
at by expanding the field operators '!/'cr(r) in terms of 
the locaHzed Wannier functions [20| . The Wannier wave- 
functions on the lowest band, wo,cr(r), are even functions, 
while the Wannier function ^^^^(r) on the p-band has a 
node in the plane normal to the coordinate axis a, giv- 
ing rise to three types of state for the spin-up fermions 
which we label x, y, and z, respectively 21]. Including 
only the leading nearest neighbor tunneling terms, we 
find the Hamiltonian for the ideal two-component Fermi 
gas in momentum space 



-ffidcal = y^ (ep.k.o- - A*<t) V'd.k.g^O.k.' 



(1) 
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where a = {Tj i}- Here the summation is over the first 
Brillouin zone, V'o u o- creates an atom in the lowest band, 
ficr are the chemical potentials |22|, and eo,k,(T and eQ.k.o- 
are the dispersions on the lowest and excited p-bands 
respectively. Furthermore, V'l k t creates a spin-up atom 
with momentum k. 

The lowest band dispersion is given by the usual ex- 
pression eo,k,<T = J2(3={x,y,z}'^'^o,<^J^i^ ~ cosk/id), where 
in principle the tunneling strength can depend on the 
spin-state as well as on direction if the lattice depth is 
different in different directions. The excited band disper- 
sions are more complex. First of all, the excited bands are 
separated from the lowest band by energy gaps Aa.sG- 
Second, since the Wannier functions on the excited band 
have a node and are antisymmetric along some axis, the 
tunneling strength for moving an atom in the direction 
orthogonal to the nodal plane is different from moving it 
in the direction along the nodal plane. The dispersions 
are then ea^k,T = ^a:BG + J2i3={x,y,z}'^-^af3 (1 -cosk/sd) , 
where Jap is the tunneling strength in the direction /3 for 
an atom which has a localized wavefunction Wo,(r). 

In ultracold gases the dominant interaction between 
unlike fermions is typically the s-wave interaction 
g J drn^(r)ni(r), where the coupling g can be expressed 
in terms of the s-wave scattering length a and atomic 
mass m as g — AttTi a/m. This interaction term can 
again be reduced into the lattice by expanding the field 
operators and keeping only the leading on-site interac- 
tion terms. This procedure gives us a contribution in 
coordinate space 

Hi = Uo ^ i'o,ui^l,i,ii'oxii'o,iA + (2) 

1 — \^x t'^'y j^z ) 

51 Y^ ^l,aV'L,TV'o,i,i^O,i,iV'a,i,T 

i— (*:c ,^y ,^2 ) CK — {3^,^,2} 

to the Hamiltonian. The coupling strength Uq between 
atoms on the lowest band is related to the scattering 
length through [/q — .g/ dr|wo,x(r)[^|wo,T(r)|^ while the 
interband couplings between the J,-atoms on the lowest 



band and the f-atoms on the excited band are given by 

In principle, the localized Wannier functions could be 
calculated numerically from the three-dimensional band 
structure, but this is unnecessary for our purposes. One 
gets reasonable analytical estimates for all the parame- 
ters of the theory by approximating the Wannier func- 
tions with harmonic oscillator states localized at each 
lattice site. General features are not sensitive to pre- 
cise numerical values of the parameters and the harmonic 
approximation gives us a handle on how various param- 
eters vary relative to one another as the lattice depth is 
changed. The analytical formulas for the parameters are 
long and not very informative and are, for that reason, 
omitted here. It is however necessary to discuss some 
general features. 

Firstly, for typical parameters the bandgaps are much 
higher than other energy scales of the problem. Second, 
in the harmonic approximation the interband coupling is 
simply Ui,a — Uq/2 and is independent of which excited 
band is involved, even for an-isotropic lattices. Third, 
on the excited bands the diagonal tunneling strengths 
Jaa have an opposite sign to the off-diagonal strengths 
(as well as tunneling strengths on the lowest band) and 
also their magnitude is much larger than those for the off- 
diagonal hopping strengths. This is a simple consequence 
of the Wannier function being an odd function of the co- 
ordinate normal to the nodal plane, as well as having a 
wider extension along that direction. This is also implied 
by the familar bandstructures in a one-dimensional sys- 
tem where ther lowest band has a minimum at fc = 
while the first excited band has a minimum at the edge 
of the Brillouin zone. 

Because of the strong anisotropy of the excited band 
tunneling strengths, the structure of the ideal gas Fermi 
surface is quite unlike that in the lowest band. In the low- 
est band, the Fermi surface is roughly spherical for low 
filling fractions, but on the excited bands Fermi surfaces 
are more sheet-like since atoms first fill the states along 
the directions perpendicular to the direction of large tun- 
neling strength. The Hamiltonian we derived is very rich 
and the possible phases that can occur in different param- 
eter regimes are numerous. We now proceed to discuss 
a few applications. First, we discuss anti-ferromagnetic 
phases with repulsive interactions and then we proceed 
to discuss the paired phases with attractive interactions. 



III. REPULSIVE INTERACTIONS 

The order parameter for an antiferromagnetic state 
with an ordering vector Q is defined as Aa^Q^a = 

V^~^ X]k('/'a k-i-Q o-V'a.k.o-), whcrc V is the dimensionless 
volume, i.e., the total number of sites in the system. In 
principle, two types of antiferromagnetic state are con- 
ceivable in the repulsive case: either one with checker- 
board symmetry, so that Q = (tt, tt, tt) in Cartesian co- 
ordinates, or a striped phase where symmetry is broken, 



Q = (tt, 0,0). However, as we shall see, the striped 
phase is found to be energetically unfavorable. Follow- 
ing Ref . [231 , the interaction part of the Hamiltonian is 
written in a mean-field approximation as 

ak 

+ Ui.arial (■0o^kiV'O,ki + -^J^k+Q.i^O.k+Qa) 
+ C^l,a-4o,Q,x (^V'l^k+Q,T^"kT +h.C.j 

- VYUl,a{noina'\ + Ao^QlAaQ^) , (3) 

a 

where the primed sum extends only over the reduced Bril- 
louin zone (RBZ), defined such that the points {k, k -|- 
Q|k £ RBZ} make up the entire first Brillouin zone. The 
quantities hqi and Uai are the densities of the compo- 
nents occupying the respective orbitals. A Bogoliubov 
transformation diagonalizes the Hamiltonian, where- 
after the free energy r2(ylo,Q,i'^2:,Q,T' ^y,QJ' ^2,Q,t) = 
— fcsTlogjTr exp[— /3(i?idoai + H{)]} can be computed. 
Here, f3 = 1/(A:bT), and T is the temperature. For defi- 
niteness, we assume that both components experience the 
same lattice potential and that the lattice potential has 
the same depth in all directions. Parameters are calcu- 
lated for ^"K atoms in a lattice of depth Vq = 8 Eji, where 
Eji is the recoil energy of the atoms when they absorb a 
photon (of wavelength 826 nm). In that case, we obtain 
Jo = 0.012£:7?, Jaa = -0.15£:fl, and J^p = 0.0028^^^ if 
a ^ (3, which demonstrates our earlier qualitative argu- 
ments about the magnitude of the tunneling terms. For 
a = 174 Bohr radii, Uq = QMEr ^. We assume that 
using a magnetic field a can be tuned to arbitrary posi- 
tive and negative values [l|. 

As a is increased, the system undergoes a first-order 
transition from the normal state, with AaQa = for 
all a, to a checkerboard state where AaQa 7^ and 
Ax,ci^ = ^y,Q,T = ^z,Q,T- The order parameter A^^q-^ is 
plotted in Fig.[TJ The magnitude of ^o,Q,i displays a sim- 
ilar behavior. The chemical potentials are in the calcula- 
tion fixed to ^x = 6Jo,x and /i| = /^x,bg + %Jxy + Jxz), 
respectively; this is in the noninteracting limit close to 
the value for half-filling. The phase diagram is found to 
be very insensitive to the exact values of the chemical po- 
tentials; in the antiferromagnetic state, the densities are 
locked to the values ni = 0.5 and n| = 1.5, respectively, 
so that for the majority component, the mean occupa- 
tion of the orbital is unity and the combined occupa- 
tion of the x, y, and z orbitals is 0.5. If the fernuons 
are confined in a magnetic trap, the state of the system 
is to a good approximation given by a local-density ap- 
proximation where the local chemical potential is given 
by the sum of the chemical potential and the negative 
of the trapping potential. The fact that the state of 
the system is insensitive to the chemical potential then 
means that without much fine-tuning, the antiferromag- 




FIG. 1: Magnitude of the antiferromagnetic order parameter 
^ajQT for the majority component in the checkerboard state, 
Q — (tTjTTjTt), as a function of temperature and scattering 
length for a gas of '"^K atoms in a lattice with strength Vb = 
8-Bh, where Er is the recoil energy. The scattering length a 
is given in units of the Bohr radius as ■ 



netic state can be made to occupy a large area in the 
trap, with corrections only toward the edges where the 
density drops to zero. 

A fully antiferromagnetic state, where the occupations 
of neighboring states in each spin state alternates be- 
tween and I, has an order parameter of magnitude ex- 
actly equal to 0.5. It is seen in Fig. [1] that the system 
in the checkerboard state is always close to this limit be- 
cause of the strong repulsive interactions. The absence 
of a striped phase in these calculations is not conclusive 
proof that such a state is always thermodynamically un- 
favorable. On the contrary, it is conceivable that striped 
phases could show up, e.g., in anisotropic lattices or in 
lower dimensions. However, for the parameter ranges 
that we investigated here, the checkerboard state always 
has the lower free energy. 



IV. ATTRACTIVE INTERACTIONS 

In order to study BCS-type paired states with attrac- 
tive interactions (i.e. negative a), we introduce aux- 
iliary (pairing) fields Aq = t/o(V'o,i,iV'o,i,T) ^'^'^ ^a = 
Ui,a{'4'o,i,i'4'a,i,'\) which we use to decouple the interac- 
tion terms in the usual way. In this way we find the 
interaction term of the mean-field Hamiltonian 

i 

+ ^ A„vI,i,TV'ai,x + A;i/;o,i,iV'a,i,T - |Aa| Vf^i,a 

a.'i 

The Hamiltonian is then diagonalized with a canonical 
transformation and the grand potential computed and 
minimized in order to find the state that is physically 
realized. 



LU 



10 
9 
8 


X 


10 


" 














































1 


^^^1 






















^^^^1 






































7 
























































































































d: 6 














































J 
































































































^^^^^^^^^^^^^^^^H 


m b 
















































^^^^^^^^^^^^^^1 




















































^^^^^^^^^^^^^H 


4 


















































^^^^^^^^^^^^^^^1 




















































^^^^^^^^^^^^H 






















































^^^^^^^^^^^^1 


3 




















































^^^^^^^^^^^^^^1 






















































^^^^^^^^^^^H 
























































^^^^^^^^^^^H 


2 






















































^^^^^^^^^^^^^1 
























































^^^^^^^^^^^H 
























































^^^^^^^^^^H 


1 






















































^^^^^^^^^^^^H 
























































^^^^^^^^^^H 
























































_^^^^^^^^^| 



1.3 



1.4 



1.5 1.6 

n^(T=0) 



1.7 



FIG. 2: Phase diagram with a fixed scattering length a = 
— 80a_B and the minority component chemical potential which 
corresponds to half-filling for the ideal system at T = 0. The 
a;-axis shows the filling fraction of the ideal gas of majority 
atoms at T = Shading: symmetric state = light, on-axis 
state = gray, normal state=dark. 



In Fig. [2] we show an example phase diagram when the 
minority component chemical potential is fixed to the 
value corresponding to half filling for an ideal Fermi gas 
at T = and the majority component chemical potential 
and temperature are varied. The lattice parameters were 
chosen the same as before. For the assumed symmetric 
lattice we can identify three different phases: the normal 
state where all pairing fields vanish, the on-axis state 
where only one A^ ^ 0, and the symmetric state where 
all pairing fields are equal, i.e., A^; = Aj^ = A^. The 
last one of these states dominates at low temperatures 
and for lower majority component filling factors. The 
on-axis state can be favorable at low temperatures and 
somewhat higher filling factors, while the normal state is 
favorable elsewhere. It should be noted, that whether or 
not the on-axis state appears depends also on the cou- 
pling strength. If the coupling is increased, the symmet- 
ric state with equal pairing fields occurs on a larger part 
of the phase diagram. 

/^From our computations (not shown in the figures) we 
observe the presence of a critical coupling strength before 
pairing can take place. The reason for the critical cou- 
pling strength has to do with a different structure of the 
Fermi surfaces of the majority and minority atoms. For 
identical atoms at the lowest band the Fermi surfaces can 
be perfectly matched for zero polarization, but this is no 



longer true when one of the components occupies states 
on the excited p-bands. In this case sufficiently strong 
coupling is required to counteract the effect of the mis- 
matched Fermi surfaces. Just above the critical coupling 
strength the on-axis state has a slightly lower energy than 
the symmetric state while for stronger couplings the lat- 
ter state is favored. 



V. CONCLUSIONS 

Optical lattice experiments usually include a parabolic 
trapping potential. When one starts to increase the num- 
ber of majority atoms while keeping the number of mi- 
nority atoms fixed, the higher bands will not be occu- 
pied straight away. At first, the cloud of majority atoms 
spreads out in the harmonic trap and their filling factor 
in the optical lattice remains less than one. However, 
roughly at distances larger than R the energy muj^R^ /2 
due to the trapping potential with frequency ojt becomes 
larger than the bandgap. When this happens, it is favor- 
able for the atoms to start filling the excited band(s) in 
the center of the trap. This implies that in a trapping en- 
vironment the phases discussed in this paper can occur in 
the center of the atomic cloud and that this center will be 
surrounded by a cloud of majority atoms occupying the 
lowest band. The presence of the pairing gaps and anti- 
ferromagnetic ordering would be observable, for example, 
through noise-correlation experiments [25|, |26|, |27| . In the 
noise correlation experiments the structure of the corre- 
lation peaks will depend on the symmetry of initial Wan- 
nier functions of the atoms prior to free expansion and 
can therefore be used to distinguish various phases. Pair- 
ing gaps could also be observed through radio-frequency 
sprectroscopy [M [H [33| . 

In this article we derived a theory for the two- 
component polarized fermions in an optical lattice when 
also the lowest excited p-bands are occupied. Based 
on this theory we studied anti-ferromagnetic phases as 
well as mean-field BCS-type theory with several possi- 
ble order parameters. We outlined the expected phase 
boundaries for anti- ferromagnetic phases, whose proper- 
ties need to be studied in greater detail in the future. 
In the attractive case we assumed BCS order parameters 
which do not break translational symmetry, thus post- 
poning the investigation of states which do break the 
translational symmetry in a lattice [7|, |3l|. Interesting 
physics is also expected when the minority component 
starts to populate the excited bands. In this case one 
can expect competition between inter-band and intra- 
band effects. 
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